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We solve the attractive Hubbard model for arbitrary interaction strengths within dynami- 
cal mean-field theory. We compute the transition temperature for superconductivity and ana- 
lyze electron pairing in the normal phase. The normal state is a Fermi liquid at weak coupling 
and a non-Fermi liquid state with a spin gap at strong coupling. Away from half-filling, the 
quasi-particle weight vanishes discontinuously at the transition between the two normal states. 
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Fermi systems with a weak attractive interaction are 
Fermi liquids which undergo a phase transition into a 
superconducting state via the condensation of weakly 
bound Cooper pairs at a low critical temperature T c . 
For a long time this weak coupling route to supercon- 
ductivity, which is well described by the BCS mean-field 
theory H, was distinguished almost dogmatically from 
Bose-Einstein condensation of bosons, although in both 
cases a phase transition characterized by the same U(l) 
symmetry breaking occurs. Leggett was the first to 
point out that BCS superconductivity can be continu- 
ously connected to Bose-Einstein condensation of tightly 
bound pairs by increasing the two-particle attraction be- 
tween the fermions from weak to strong coupling. The 
size of the Cooper pairs shrinks continuously until spa- 
tially well separated bosons form, which undergo Bose 
condensation at a sufficiently low temperature. Nozieres 
and Schmitt-Rink [|| have extended Leggett 's analysis to 
lattice electrons and finite temperatures. Based on phys- 
ical insight gained from a discussion of the weak and 
strong coupling limits and an approximate (T-matrix) 
treatment of the intermediate regime they concluded that 
the evolution from weak coupling to strong coupling su- 
perconductivity is indeed smooth. 

The interest in the intermediate regime between the 
BCS and Bose-Einstein limits increased considerably af- 
ter the discovery of high-temperature superconductors, 
which are characterized by Cooper pairs whose size is 
only slightly bigger than the average electron distance 
j|. Much recent work has therefore been dedicated to 
the theory of Fermi systems with attractive interactions 
of arbitrary strength It was shown that sufficiently 
strong attraction or reduced dimensionality can lead to 
energy gaps even in the normal phase, which have been 
related to pseudogap phenomena in the cuprate super- 
conductors ||. 

In this work we analyze the formation of pairs in a 
Fermi system with arbitrary attractive interactions by 
solving the attractive Hubbard model within dynamical 
mean- field theory (DMFT) . We show that the normal 
state is a Fermi liquid at weak coupling and a non-Fermi 
liquid state characterized by bound pairs and a spin gap 



at strong coupling, in qualitative agreement with Quan- 
tum Monte Carlo (QMC) studies of the two- and three- 
dimensional attractive Hubbard model |§Hnfl- At very 
low temperatures the transition between the Fermi liq- 
uid and the normal paired state is discontinuous. 

In standard notation the Hubbard model for lattice 
fermions with a nearest neighbor hopping amplitude —t 
and a local interaction U is given by 

ff = - f EE c t c > + ^E^t^i ■ (i) 

(U> " 3 

The attractive (U < 0) Hubbard model is a supercon- 
ductor below a certain critical temperature T c (U,n) > 
for all U at any average density n, if the lattice dimen- 
sionality is above two At half- filling (n = 1) the 
usual J7(l) gauge symmetry becomes a subgroup of a 
larger 5*0(3) symmetry, and the superconducting order 
parameter mixes with charge density order. In two di- 
mensions one expects a Kosterlitz-Thouless phase at low 
temperatures for all U < and n ^ 1 |p"lj| . 

In the weak coupling limit C/ — > and dimensions d > 2 
the attractive Hubbard model can be treated by BCS 
mean- field theory p|JTl|]. In the strong coupling limit 
U — > — oo the low energy sector of the model (excitation 
energies <C |?7|) can be mapped onto an effective model 
of hard core lattice bosons with a hopping amplitude of 
order t 2 /U and a repulsive nearest neighbor interaction 
of the same order . These bosons undergo Bose con- 
densation in d > 2 dimensions and a Kosterlitz-Thouless 
transition in two dimensions (for n ^= 1) at a critical 
temperature of order t 2 /\U\. 

For nearest neighbor hopping on a bipartite lattice the 
particle-hole transformation of spin-"f fermions 

where r/j — 1 (—1) for j on the A-sublattice (B- 
sublattice), maps the attractive Hubbard model at den- 
sity n onto a repulsive Hubbard model at half-filling with 
a finite average magnetization m = 1 — n |TT| . We will 
use this relation to compare with results known for the 
repulsive Hubbard model. 
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We have solved the attractive Hubbard model within 
DMFT (7). In contrast to other (simpler) mean-field ap- 
proaches, DMFT provides an exact solution of the model 
in the limit of infinite lattice dimensionality | fl2f , since it 
captures local fluctuations exactly. 

To convince ourselves that DMFT is a suitable ap- 
proach for the weak to strong coupling crossover problem, 
let us first consider the limits. At weak coupling DMFT 
captures the complete BCS physics, since it contains the 
Feynman diagrams contributing to the BCS mean-field 
theory. At strong coupling, where the attractive Hub- 
bard model maps to the hard core Bose gas, DMFT re- 
duces to the standard mean-field theory of the hard core 
Bose gas Hence, Bose-Einstein condensation of pre- 
formed pairs is obtained at a critical temperature of order 
t 2 /\U\ at large \U\. 

Within DMFT the fluctuating environment of any lat- 
tice site is replaced by a local but dynamical effective field 
Go @. The mean-field equations involve the calculation 
of the propagator G(t) = — {Tc a {T)c\, (0)} of an effec- 
tive single-site Hubbard model coupled to the dynamical 
field Qq, and a self-consistency condition relating G to 
the local propagator of the full lattice system. The lat- 
tice structure enters only via the bare density of states 
(DOS), as long as the translation invariance of the lat- 
tice is not broken. We have used the particularly simple 
self-consistency equations 

Qq 1 ^) = iu + (i- (e /2) 2 G(iuj) (3) 

valid for a half-ellipse shaped density of states Dq(e) — 
— e 2 . Any other bounded DOS would yield qual- 
itatively similar results. In the following we will set 
eo/2 = 1. Susceptibilities such as the pairing and the spin 
susceptibility can also be computed from expectation val- 
ues of operator products within the effective single-site 
problem. The DMFT equations can also easily be ex- 
tended to superconducting or other symmetry broken 
phases 0. In this work we focus however on normal 
state properties. 

The effective single-site problem cannot be solved an- 
alytically. We have solved it numerically by discretizing 
the imaginary time interval and computing expectation 
values via the standard Hirsch-Fye algorithm |l4|] . 

We now present and discuss results from our DMFT 
calculation at quarter-filling (n = 1/2). We do not ex- 
pect that the results depend qualitatively on the density 
in the attractive Hubbard model, as long as n is finite. 
Only the particle-hole symmetric case n = 1 is special 
due to its larger symmetry group. Quarter-filling is well 
below half-filling but still high enough to see collective 
many-body effects, which are not obtained in the low- 
density limit. 

The critical temperature T c for the onset of supercon- 
ductivity shown in Fig. 1 has been obtained from the 
pairing susceptibility in the normal phase, which diverges 
as the temperature approaches T c from above. The 
critical temperature at half-filling, where the supercon- 



ducting state is degenerate with a charge-density wave 
state, has been computed within DMFT already earlier 
Jl5| . Correlations of course suppress T c with respect to 
the BCS result. At strong coupling T C (U) approaches 
the 1/U behavior described by the hard core Bose gas 
limit. At intermediate coupling strengths, T C (U) varies 
smoothly, as predicted already by Nozieres and Schmitt- 
Rink A comparison with a recent result for T c ob- 
tained by combining the DMFT with a self-consistent T- 
matrix approximation (TMA) |ll| shows that the latter 
approximation reproduces the correct qualitative behav- 
ior of T C (U), but fails quantitatively. 
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FIG. 1. Critical temperature T c as a function of the cou- 
pling strength \U\ at quarter-filling within DMFT, compared 
to T c obtained from BCS theory and from the T-matrix ap- 
proximation (TMA), respectively. 

In the following we discuss the weak to strong coupling 
crossover in the normal phase. We ignore the supercon- 
ducting instability and study normal state solutions of 
the DMFT equations also below T c . Of course these so- 
lutions do not minimize the free energy, but they could 
be stabilized by the field energy of a sufficiently strong 
external magnetic field. 

At weak coupling the normal state of the system is 
a Fermi liquid with fermionic quasi-particle excitations. 
Besides numerical evidence this follows fl7|| from the an- 
alyticity of weak coupling perturbation theory for the 
effective single-site problem. By contrast, at sufficiently 
strong coupling \U\ 3> eo and zero temperature all par- 
ticles arc bound in pairs, because a small kinetic energy 
cannot overcome a finite binding energy. Only short- 
ranged virtual breaking of local pairs occurs. At low 
finite temperatures T <C \U\ only an exponentially small 
fraction of pairs dissociates. 

A good measure for local pair formation is the density 
of doubly occupied sites rid = ( n jT ■ For an uncorre- 
cted state the density of doubly occupied sites is simply 
the product of the average density of up and down spin 
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particles, i.e. n° d = n± = (n/2) 2 . An attractive inter- 
action enhances rid- In the limit of infinite attraction all 
particles are bound as local pairs such that rid — > n/2. 
In Fig. 2 we show results for rid(T) for various U. 
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FIG. 2. Density of local pairs as a function of tempera- 
ture for various coupling strengths U at quarter-filling. 

For T — > oo the density of doubly occupied sites tends 
to (n/2) 2 , corresponding to an uncorrelated state. For 
decreasing temperature tid(T) first increases as a conse- 
quence of the attractive interaction. For small or mod- 
erate U, however, rid(T) slightly decreases again at low 
temperatures. This effect, which has also been obtained 
in the combined DMFT + TMA calculation @ , can be 
attributed to the kinetic energy, which tends to dissoci- 
ate pairs if the attraction is not too strong. Note that in 
the pairing regime for stronger U the upturn in nd(T) at 
low temperatures is missing. The kinetic energy is not 
able to unbind pairs any more. A completely analogous 
(particle-hole transformed) behavior has been found in 
the DMFT solution of the repulsive Hubbard model at 
half-filling §. 

The absence of fermionic quasi-particles in the pair- 
ing state at strong coupling also leads to a pronounced 
spin gap. In Fig. 3 we show our DMFT results for the 
temperature dependence of the spin susceptibility Xs, for 
various coupling strengths. For a weak attraction the 
spin susceptibility increases monotonously for lower tem- 
peratures and then saturates at a finite value for T — > 0, 
as expected for a Fermi liquid. For strong coupling, how- 
ever, Xs decreases rapidly at low temperatures, as ex- 
pected for a system where spin excitations are gapped. 
This gap, which has also been seen in QMC simulations 
of the two-dimensional [||| and three-dimensional |lC| ] 
Hubbard model, is clearly due to the binding energy of 
pairs in the non- Fermi liquid state forming at strong cou- 
pling. For a moderate attraction, pseudogap behavior 



seems to set in at intermediate temperatures, but for 
T — > the susceptibility tends to a reduced but non-zero 
value. We expect that this behavior reflects the presence 
of a narrow quasi-particle band in the system, equivalent 
to the one known for the repulsive model near the Mott 
transition . The selfconsistent TMA fails to yield spin 
gap behavior in high dimensions fll6| |, and yields only a 
rather weak suppression of Xs in two dimensions [18| . 
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FIG. 3. Spin susceptibility Xs as a function of temperature 
for various coupling strengths U at quarter-filling. 

Since the Fermi liquid state at weak coupling is qual- 
itatively different from the bound pair state at strong 
coupling, there must be a sharply defined pairing transi- 
tion at some critical attraction U c at least in the ground 
state. At quarter-filling we can estimate from our data 
C7 c «-2.5e . 

To see how the Fermi liquid breaks down upon increas- 
ing the attraction strength, we have computed the renor- 
malization factor Z(T) = [1 — S(iwo) / Wo] -1 , where £ 
is the self-energy and too — irT the smallest (positive) 
Matsubara frequency at temperature T. In Fig. 4 we 
plot Z(T) as a function of T for various U at quarter- 
filling (left) and at half- filling (right). At quarter- filling 
Z(T) extrapolates to a finite positive value Z in the limit 
T — > 0, for any U . In the Fermi liquid phase Z has physi- 
cal meaning, being the spectral weight for quasi-particles, 
the Fermi edge discontinuity in the momentum distribu- 
tion function and, within DMFT, also the inverse mass 
renormalization. This meaning is of course lost in the 
bound pair state. Note that the finiteness of Z for all 
U does not imply that the system is a Fermi liquid for 
arbitrary interactions! A simple calculation shows that 
Z is finite even in the atomic limit t = 0, where the sys- 
tem is obviously not a Fermi liquid. An exception is the 
half- filled case, where Z — > for U — > U c , and Z = 
for all \U\ > \U C \ (see Fig. 4) and in the atomic limit. 
For half-filling the continuous vanishing of Z has been 
analyzed in much detail for the Mott transition in the 
repulsive Hubbard model M, which is equivalent to the 
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pairing transition in the attractive model by virtue of the 
particle- hole symmetry (2). We thus conclude that the 
Fermi liquid phase disappears with a finite Z at electron 
densities n ^ 1, i.e. the quasi-particle weight disappears 
discontinuously at the pairing transition. 
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FIG. 4. Z(T) as a function of temperature for various U at 
quarter-filling (left) and half-filling (right). 

It is instructive to consider the low density limit n — > 
for comparison. In that limit the bound pair state is 
stable once the attraction exceeds the threshold for two- 
particle binding [/?. For \U\ < \U°\ no bound states 
exist, and the particles move essentially freely, due to 
the low density, and Z is almost one even close to the 
transition. 

Using the Hirsch-Fye algorithm it is hard to determine 
unambiguously whether there is a sharp phase transition 
at some critical U C (T) also at sufficiently low finite tem- 
perature or only a smooth (though steep) crossover, since 
the computation time increases rapidly at low tempera- 
tures. However, one can find an answer to this question 
by exploiting the equivalence of the attractive Hubbard 
model at generic densities and the half-filled repulsive 
Hubbard model with a finite magnetization. The latter 
model has been analysed earlier within DMFT by Laloux 
et al. [^9|, who solved the DMFT equations with an ex- 
act diagonalization algorithm which is more efficient than 
the Hirsch-Fye algorithm at low temperatures. The re- 
sults of their work imply that at very low temperatures 
a first order transition occurs in the attractive Hubbard 
model between a thermally excited Fermi liquid state and 
a thermally excited bound pair state. 

In summary, we have shown that DMFT yields a tran- 
sition from a Fermi liquid state at weak coupling to a 
non- Fermi liquid state with a spin-gap at strong coupling 
in the attractive Hubbard model. Spin-gap behavior for 
strong attraction is obviously governed by short-range 
pair correlations which are captured by DMFT. Long- 
range superconducting fluctuations characteristic for low- 
dimensional systems are crucial for normal state gaps 
only at weaker coupling. 
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